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Homotopy Analysis to Heat and Mass Transfer of 
MHD Flow of a Viscous Fluid over a Moving 
Vertical Plate in a Porous Medium with Viscous 

Dissipation 
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Abstract — The present paper deals with the MHD boundary 
layer flow over a linearly moving porous semi infinite vertical 
plate taking suction and viscous dissipation into account. The 
fluid considered is viscous, incompressible and electrically 
conducting. After transferring the governing equations into 
ordinary differential equations using suitable dimensionless 
variables, analytical solutions are generated using homotopy 
analysis method (HAM). The effects of various parameters on 
dimensionless velocity, temperature and concentration profiles 
are presented in the form of graphs and tables. HAM results are 
in good agreement with the results available in the literature. 


Index Terms — Porous Media, Viscous Dissipation, MHD 
Flow, Moving Vertical Plate, Homotopy Analysis 
Method(HAM). 


I. INTRODUCTION 

In recent years, there are number of studies on convection 
flow and heat transfer in saturated porous media due to their 
wide ranging applications in engineering field like heat 
exchanger devices, geothermal and geophysical engineering, 
petroleum reservoirs, underground disposal of nuclear waste 
and others. The vertical free convection boundary layer flow 
in porous medium owing to combined heat and mass transfer 
has been studied by Bejan and Khair [1], Lai and Kulacki [2] 
studied the coupled heat and mass transfer by natural 
convection from vertical surface in porous medium. Kim and 
Vafai [3] have analyzed the buoyancy driven flow about a 
vertical plate for constant wall temperature and heat flux. 

Helmy [4] worked on MHD unsteady free convection flow 
past a vertical plate embedded in a porous medium. Raptis et 
al. [5] constructed similarity solutions for boundary layer near 
a vertical surface in porous medium with constant temperature 
and concentration. Thermal radiation and buoyancy effects on 
hydromagnetic flow over an accelerating permeable surface 
with heat source or sink was investigated by Chamka [6], Yih 
[7] observed free convection effect on MHD coupled heat and 
mass transfer of a moving permeable vertical surface. 
Soundalgekar [8] has discussed the viscous dissipation effect 
on unsteady free convection flow past an infinite vertical 
porous plate with constant suction. 

In this study, HAM [10, 11] is adopted to find the analytical 
solutions for velocity, temperature and concentration profiles. 
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Convergence of the obtained solutions is explicitly shown. 
Velocity, temperature and concentration profiles are drawn 
for various values of flow parameters like magnetic 
parameter M , buoyancy parameters Gr and Gc , suction 

parameter F , permeability parameter K , Schmidt 
number Sc and Eckert number Ec . The results are 
discussed in detailed, compared with the available results in 
literature and are in good agreement. 

Governing equations 

Consider the free convection effects on MHD boundary 
layer flow of a viscous incompressible fluid over a linearly 
started porous vertical semi infinite plate embedded in a 
porous medium with suction and viscous dissipation. The 
X -coordinate is taken along the plate in ascending direction. 
The y -coordinate is taken normal to the plate. The velocity 
of the fluid far away from the plate surface is assumed to be 
zero for a quiescent state fluid. The variations of surface 
temperature and concentration are linear. A uniform magnetic 
field is applied in the direction perpendicular to the plate. The 
transverse applied magnetic field and magnetic Reynolds 
number are assumed to be very small, so that the induced 
magnetic field. Hall effects and Joule heating are negligible. 
Under these assumptions, along with the Boussinesq 
approximations, the boundary layer equations describing this 
flow are 


du dv 
dx dy 


( 1 ) 


du du 
u — tv — 
dx dy 



oB 


+80t {t-tJ+sMc-cJ- 



(2) 


dT dT d 2 T v ( du '] 2 

u hv — = a r-H — 

dx dy dy~ pc p {dy) 


dC dC n 
u — + v — = D„ 
dx dy 


d 2 C 

dy 2 ' 


The corresponding boundary conditions are 


(3) 

(4) 


u = Bx, v = V, T =T = T +ax, C = C w =C x +bx at y = 0, 

' ' W 00 ' W 00 J > (5) 

u-> 0, T ->T X , C — ) C x asy^tw, 

where U and V are the velocity components along the X 
and y axes respectively. T and C are the temperature and 
concentration in the boundary layer, a and b denote the 
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stratification rate of the gradient of ambient temperature and 
concentration, V is the kinematics viscosity, p is the 

density, G is the electric conductivity of the fluid, p T is the 

volumetric coefficient of thermal expansion, /? is the 

volumetric coefficient of concentration expansion, C x is the 

free stream concentration, B 0 is the magnetic induction, 

D m is the mass diffusivity, g is the acceleration due to 

gravity, K ' is the permeability of the porous medium and 

OC is the thermal diffusivity, c p is the specific heat at 

constant pressure. 

Now using the non-dimensionless parameters and variables 
reported in [8] 



V(x, v) = xjvB f(rj). 


m= 


T~T X 

T V ~T X 


, <Kn)= 


c-c x 

c v -c x 


( 6 ) 


where l//(x,y) is the stream function, f (p) is a 
dimensionless stream function, 0{rj) is a dimensionless 
temperature of the fluid in the boundary layer region, 
< j>(j ] ) is a dimensionless species concentration of the fluid in 
the boundary layer region and I] is the similarity variable. 
Taking the Cauchy -Riemann equations 




(7) 


the conservation of mass, equation (1) is identically satisfied. 
Substituting equation (7) in equations (2) to (5), we obtain 

/"'+/ f'+Gr9 + Gc <f>-K f'=0, (8) 

0"+Pr f 0'-Pr f ' 0 + Ec f" 2 = 0, (9) 

</>"+Scf0'-Scf'0 = 0. (10) 

The boundary conditions in non-dimensional form are 

f = ~F w , /'=!, 0 = i <t> = 1 at 17 = 0, 

/ '= 0, 9 = 0, <f> = 0, as i] -» oo, 


where 


M = 


pB 

gfi T (T w -TJ . 


the magnetic 


parameter, Gr = w . — is the local temperature 

xB 2 

Grashof number, Gr = — ^ ‘ — is the local 

xB 2 


concentration Grashof number, K = is the 

K'xB 

v 

Permeability parameter, Pr = — is the Prandtl 

a 

„ B 2 x 2 

number, Ec = is the Eckert 

c p (T w -TJ 

v V 

number. Sc = is the Schmidt number, F, = , is 

the suction parameter. 

Solution using HAM 

To solve the dimensionless equations (8) to (10) together 
with the boundary conditions (11) analytically using HAM, 
the initial approximations and auxiliary linear operators are 
chosen as 

0M = e-\ 

^0(7) = ^, 

L 2 ( 0 ) = 6"-0, 

which satisfy 

L l (c i +C 2 e’ 1 +C 3 g~’)=0, 

L 2 (c 4 e n +C 5 e^)= 0, 

L 3 (c 6 e n + C 1 e ' 1 ) = 0, 

where C. (/ = 1 to 7) are the arbitrary constants. 

If p e [0, 1] is the embedding parameter, and h 3 are 

the non-zero auxiliary parameters, the following equations are 
constructed. 


Zeroth order deformation equations 

(1 - p) U (/ ( 7 ; p) - fo ( 0 )) = p ft 1 Ni [f fo; p\ p\ ^ p)\ ( 1 2 ) 


(l - p)L 2 (0(p; p) - 0 Q (tj)) =ph 2 N 2 [ f(rj ; p), 0{rp, p)\, (1 3) 

(l-p)L 3 (</>{wp)-</>o{n))=P03N 3 [f(p-,p),0(r/;p)l (14) 


f(0-p) = -F w , 

r( 0 ;p)=l /'(°o;p) = 0, 


0(o-,p)=i, 

6>(oo; p) = 0, 

(15) 

^(0; p) = 1, 

(p{ 00; p) = 0, 
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N i [f(p; p\ 0{rj\ p), (f>{jp, p )] = 


d 3 f(y,p) 

dp 3 


+ f(wp) 1 - 


d 2 ffap) (Sf{rr,p)\ _ M df{p\p) 


dp- 


dp 


dp 


+ Gr0(r)\ p)+ Gc</)(rp, p)- K 


Sf{p\ p) 
dp 


dp “ dp dp y dp~ 

[ f(p ; p), (/){ n ; p)]= 9 + Sc fipr, p)M^i _ 5c g /(^p) 

dp- dp dp 

Taking p — Oand p — 1 , we have 

f{r/,0)=fM /(7;l) = /tt 

0(p-O)=0 o (p), 0(p-l)=0{p), 

</>(n-,o) = <I> q (p), </>{p-,i) = </>{p), 

when /? increases from 0 to 1 then / (//; p^,0(r/', /?)and <^(//; p) vary from initial approximations to the final solutions. 

Expanding f{p', p),0{p', p)and(f>(p; p) in Taylor’s series w.r.to p , we have 

00 

fi r l) = fo( r l) + T J fMp m ’ 


m = 1 


where 


0{p)=0 0 {p)+Y^O m {p)p m , 

m = 1 

oo 

</>(p)=<i>o( r i)+Tj<i> m ( r i)p m ’ 

m = 1 

1 0'/(i7;p) 


(16) 

(17) 

(18) 


fm{ T l) = 

0 m (p) = 


ml dp m 
1 d m 0{p\p) 


p = o 


m! dp m 

_ 1 d'Vfop) 

m! dp m 


p = o 


p=0. 


If the initial approximations, auxiliary linear operators and non-zero auxiliary parameters are chosen in such a way that the 
series (16) to (18) are convergent at p — 1, then 


co 

f(p)=fo( f i)+'ZfM 

m = 1 
oo 

0 ( T j) =0 o( T ?)+^M 


m = 1 


0(7) = $>(7)+ IX W 

m = 1 

m th order deformation equations 

Differentiating equations (12) to (14) m times w.r.to p then setting p = 0 and finally dividing with m ! , we get. 


A (/» ( 7 ) - z,„ f m—\ (n)) = * 1 ^ ( 7 ), 

(19) 

L 2 (<?„, ( 7 ) - ( 7 )) = * 2 K ( 7 ), 

(20) 

L 3 ( 7 ) - Z,„ &_1 ( 7 )) = h 3 ( 7 ) 

(21) 


with the following boundary conditions 
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where 


/.( 0 )= 0 , /„( 0 )= 0 , /»= 0 , 
A(o)=o, «,„(») =o. 


(22) 


m — 1 


m - 1 


R L (>l) = / m -l + fi -Hfm-l-i fi ~ M f'm- 1 + ^ ^ + Gc^_, - K , 

i = 1 i=l 

m-l m— 1 m-l 

C, + PrSx,-.-, 


1=1 


i = 1 


i=l 


m-l m— 1 

^ l 7 /) = <Cl +5C X ~ SC Yj f m-X-i </>’ 

i = 1 i=l 

[0, 777 < 1, 

Xh = I , , 

[1, 777 >1. 

Nowtofind f m {n\ 0 m {rj) and (f> m (//) for 777 > 1, solve the equations (19) to (21) using MATHEMATIC A by employing 
the boundary conditions (22). 


Convergence of HAM 

As pointed by Liao [12], the convergence of solution 
series depends upon the choice of the initial approximations, 
the auxiliary linear operators and on the non-zero auxiliary 
parameters. Once if the initial approximations and the linear 
operators have been selected, the convergence of the solution 
series will depend upon the non-zero auxiliary parameter. 
Proper values of the parameters fl x ,tl 2 Slid tl 3 can be found 
by so-called tl -curves. According to Fig. 1, the convergent 
region of /"( 0), #'(0) and 0'(O) is [-2.0, 0.0]. In this 
paper we choose fl x =tl 2 =fl i = —1.0. Convergence of 

f ' ' (0), 0 ' (O) and (J ) ' (o) for different orders of 
approximation is given in Table 1. This table shows that the 
convergence is achieved at 30 th order approximation. 

II. Results and discussions 

The effects of various parameters such as magnetic 
parameter M , buoyancy parameters Gr and Gc , suction 

parameter F , permeability parameter K , Schmidt 
number Sc and Eckert number Ec on velocity, temperature 
and concentration fields are shown graphically in Figs. 2-22. 
In this study Prandtl number Pr is considered to be 0.72 
(air), Schmidt number Sc — 0.24,0.68,0.78, 2.72 . To 

ensure the accuracy of the present results, comparison is made 
with the existing results in the literature and is given in Table 
2 . 

Figs. 2 to 4 show the effect of magnetic parameter on 
velocity, temperature and concentration profiles. Introduction 
of transverse magnetic field to an electrically conducting fluid 
develops a drag due to Lorentz force which tends to resist the 
fluid flow and thus reducing its velocity and to increase its 
temperature and concentration. 

The effects of Gr and Gc on velocity, temperature 
and concentration fields are illustrated in Figs 5 to 10. As 
shown, temperature and the concentration are decreasing with 


increase in, Gr and Gc but the velocity increases as 
Gr and Gc increase. 

Figs. 11 to 13 depict the effect of suction parameter. 
As the suction parameter F increases f \0 and (j) are 
also increase. 

Figs. 14 to 16 show the dimensionless velocity, 
temperature and concentration profiles for different values of 
permeability parameter K . It can be seen that the velocity 
profiles decrease with the increase of permeability 
parameter K and also noticed that the temperature and 
concentration profiles increase with the increase of 
permeability parameter K . 

Figs. 17 to 19 show the dimensionless velocity, 
temperature and concentration profiles for different values of 
Schmidt number Sc . It can be seen that the velocity profiles 
decrease with the increase of Schmidt number Sc . Further, it 
is observed that the temperature monotonically increases with 
the increase of Schmidt number Sc . It is seen that the 
concentration decreases as Schmidt number Sc increases. 

Figs. 20 to 22 depict the behaviour of velocity, temperature 
and concentration profiles for different values of Eckert 
number. From the figures it is clear that there is a slight 
increase in velocity and a slight decrease in concentration 
profiles with an increase in Ec . It is also observed that 
temperature increases with Ec . 

III. Conclusions 

Analytical study has been performed to observe the 
physical behavior of the velocity, temperature and 
concentration profiles. Results are presented graphically and 
analyzed. By comparing the present results with previous 
work, it is found that there in a good agreement. 
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Fig. 1: fi -curves for f (0),f?'(0)and^ (0) a t20 th 
order approximations when 

M = 0.1, Gr = 0.1, Gc = 0.1, £ = 0.1, F w = 0.1, Ec = 0.1,5c = 0.62, Pr = 0.72. 



Fig. 2: Variation of M on f ' when 


Gr = 0. 1, Gc =0.1, K = 0. 1, F w = 0. 1, Ec = 0. 1, 5c = 0.62 , Pr = 0.72. 



Fig. 4: Variation of M on (j) when 


Gr = 0. 1, Gc = 0. 1, K = 0. 1, F w = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr = 0.72. 



Fig. 5: Variation of Gr on f ' when 


M =0.1,Gc = 0.1,X' = 0.1,F w = 0.1,£c = 0.1,5c = 0.62,Pr = 0.72. 



Fig. 3: Variation of M on 0 when 

Gr = 0. 1, Gc = 0.1 K = 0. 1, F w = 0. 1, Ec = 0. 1, Sc = 0.62, Pr = 0.72. 



Fig. 6: Variation of Gr on 6 

when 


M = 0. 1, Gc = 0. \,K = 0. 1, F w = 0. 1, Ec = 0. 1, 5c = 0.62, Pr = 0.72. 
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Fig. 7: 


Variation of Gr on (f) when 


M = 0. 1 , Gc = 0. 1, £ = 0 . 1 , F w = 0 . 1, Ec = 0. 1 , Sc = 0.62 , Pr = 0 . 72 . 



Fig. 10: Variation of Gc on (f) when 

M=()\Gr=().\,K = 0. 1, F w = 0. 1 , Ec = 0. 1 , Sc = 0.62 , Pr = 0.72. 



Fig. 8: Variation of Gc on f ' when 

M = 0. 1, Gr = 0. \,K = 0. 1, F w = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr = 0.72. 



Fig. 11: Variation of F w on /' when 

M = 0. 1, Gr = 0. 1, Gc=().\,K = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr - 0.72. 



Fig. 9: Variation of Gc on 9 when 

M = 0. 1, Gr = 0. \,K = 0. 1, F = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr = 0.72. 



Fig. 12: Variation of F w on 9 when 

M = 0. 1, Gr = 0. 1, Gc = 0. \,K = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr - 0.72. 
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Fig. 13: Variation of F on (f) when 

M = 0.1, Gr = 0.1, Gc=0.U=0.l,£c = 0.1, Sc = 0.62, Pr - 0.72. 



Fig. 14: Variation of K on f ' when 

M = 0. 1, Gr = 0. 1, Gc = 0. 1, F w = 0. 1, Ec = 0. {Sc = 0.62 , Pr = 0.72. 



Fig. 15: Variation of K on 6 when 



Fig. 16: Variation of K on (f) when 

M = 0.1, Gr = 0.1, Gc = 0.1, F w - 0.1, & - 0.1, Sc = 0.62 , Pr = 0.72. 



Fig. 17: Variation of Sc on /' when 

M = 0.1, Gr = 0.1, Gc =0.1, = 0.1, ^ = 0.1, £c = 0.1, Pr = 0.72. 



Fig. 18: Variation of Sc on 9 when 


M = 0. 1, Gr = 0. 1, Gc - 0. 1, F w = 0. 1, Ec = 0. 1, Sc = 0.62 , Pr = 0M.= 0. 1, Gr - 0. 1, Gc = 0. 1, F w = 0. \,K - 0. 1, Ec - 0. 1 , Pr = 0.72. 
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Fig. 19: Variation of Sc on (j) when 

M =0.1, Gr= 0.1, Gc = 0.1, F w = 0.1,£ = 0.1, Ec = 0.1, Pr = 0.72. 



Fig. 20: Variation of Ec on /' when 

M = 0. 1, Gr = 0. 1, Gc = 0. 1, F w =ti.l,K = 0. 1, Sc = 0.62 , Pr = 0.72. 




0 2 4 6 8 10 


0 

Fig. 22: Variation of Ec on (j) when 

M = 0. 1, Gr = 0. 1, Gc = 0. 1, F w = 0. U = 0. 1, Sc = 0.62 , Pr = 0.72. 


Table 1: Convergence of the HAM solutions for different 
order of approximations when 


M = 0.1, Gr - 0.1, Gc = 0.2, F w = 0.1, Sc = 0.78, K = 0.1, Pr = 0.72, Ec = 0.1. 


Order 

-/"( o) 

-0'(o) 

-f(o) 

5 

0.940911 

0.750398 

0.82345 

10 

0.940177 

0.74813 

0.822822 

15 

0.940134 

0.747914 

0.822905 

20 

.94013 

0.747877 

0.822931 

25 

.94013 

0.747868 

0.822937 

30 

0.940129 

0.747866 

0.822938 

35 

0.940129 

0.747866 

0.822939 

40 

0.940129 

0.747866 

0.822939 


0 t c ‘ T 

0 2 4 6 8 10 

n 

Fig. 21: Variation of Ec on 9 when 

M = 0 . 1 , Gr = 0 . 1 , Gc = 0 . 1 , F w = 0. \,K- = 0. 1 , Sc = 0.62 , Pr = 0.72. 
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Table 2: Comparison of the results of — f " (0), — 0' (0) and — </>' (0) with the existing results of Ibrahim and Makinde [10] 
when Pr = 0.72, K = Ec = 0.0. 


Gr 

Gc 

M 

F 

Sc 

Ibrahim [10] 

HAM Results 



-/"( 0) 

-d\0) 


-/"( 0) 

-0'(O) 


0.1 

0.1 

0.1 

0.1 

0.62 

0.888971 

0.7965511 

0.7253292 

0.888967 

0.7965538 

0.7253252 

0.5 

0.1 

0.1 

0.1 

0.62 

0.695974 

0.8379008 

0.7658018 

0.695973 

0.8379015 

0.7658006 

1.0 

0.1 

0.1 

0.1 

0.62 

0.475058 

0.8752835 

0.8020042 

0.475071 

0.8753137 

0.8020275 

0.1 

0.5 

0.1 

0.1 

0.62 

0.686927 

0.8421370 

0.7701717 

0.686926 

0.8421389 

0.7701723 

0.1 

1.0 

0.1 

0.1 

0.62 

0.457723 

0.8818619 

0.8087332 

0.457704 

0.8818575 

0.8087296 

0.1 

0.1 

1.0 

0.1 

0.62 

1.264488 

0.7089150 

0.6400051 

1.264514 

0.7088625 

0.6400368 

0.1 

0.1 

0.1 

1.0 

0.62 

0.570663 

0.5601256 

0.5271504 

0.570658 

0.5601314 

0.5271426 

0.1 

0.1 

0.1 

0.1 

0.78 

0.893454 

0.7936791 

0.8339779 

0.893453 

0.7936779 

0.8339784 

0.1 

0.1 

0.1 

0.1 

2.62 

0.912307 

0.7847840 

1.6504511 

0.928472 

0.7931287 

1.9265042 
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